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Abstract 

Highly supercritical accretion discs are probable sources of dense optically thick 
axisymmetric winds. We introduce a new approach based on diffusion approximation 
radiative transfer in a funnel geometry and obtain an analytical solution for the 
energy density distribution inside the wind assuming that all the mass, momentum 
and energy are injected well inside the spherization radius. This allows to derive 
the spectrum of emergent emission for various inclination angles. We show that self- 
irradiation effects play an important role altering the temperature of the outcoming 
radiation by about 20% and the apparent X-ray luminosity by a factor of 2 — 3. The 
model has been successfully applied to two ULXs. The basic properties of the high 
ionization HH-regions found around some ULXs are also easily reproduced in our 
assumptions. 
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Introduction 



Processes of gas accretion onto compact objects are studied since 60s when first X-ray 
sources were discovered (Frank et al. 2002). Among several works describing the details of 
disc accretion in binary systems Shakura & Sunyaev (1973) was the most successful!. Their 
"standard disc" is still widely used to describe the thermal component of X-ray spectra of X-ray 
binaries. Standard disc model was worked out in several considerations such as high optical 
thickness, low geometrical thickness of the disc, etc. Among these assumptions one was that 
the power released in the accretion process does not exceed the Eddington luminosity. This 
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case is usually called subcritical accretion. 

However, considering various phenomena such as growth of supermassive black holes 
and nova outbursts leads to the problem of supercritical accretion. It has been extensively 
investigated since 1980s when Abramowicz proposed it as a power source for active galactic 
nuclei (Abramowicz et al. 1980). Presently super-Eddington accretion is often applied to explain 
the observational properties of Ultraluminous X-ray Sources (ULXs, see Roberts (2007) for 
observational review) that are believed to be high-mass X-ray binaries with black holes accreting 
on thermal time scale. 

The first, outflow-dominated model of supercritical accretion, has been developed by 
Shakura and Sunyaev in the paper mentioned above. Authors assumed that all the inflowing 
gas above the critical accretion rate is being ejected from the disc in the form of a wind. 
Another version of supercritical regime is based on the relaxation of the locality condition 
in "standard model" which leads to advective slim discs (Abramowicz et al. 1980) or Polish 
doughnuts (Abramowicz et al. 1978; Kozlowski et al. 1978; Abramowicz 2004). In reality, both 
processes work simultaneously. 

A comprehensive model taking into account both for advection and ejection was recently 
developed by Poutanen et al. (2007). Authors consider the structure of the disc in radial 
direction and estimate three characteristic temperatures relevant for the outcoming radiation 
(inner disc temperature, temperature at spherization radius and effective temperature of the 
wind photosphere) but do not study the processes of radiation transfer in the wind and do not 
calculate the outcoming spectra. 

Obviously, at high accretion rates the observational properties of the accretion flow are 
governed mostly by radiative transfer in the outflowing wind. The observational appearance 
of the pseudo-photosphere of the optically-thick wind of a supercritical accretion disc was 
considered by Nishiyama et al. (2007). The authors however assumed the optically-thick part 
of the wind fully adiabatic not taking into account radiative energy transfer that is expected 
to be important in the flow (see section 2.1). 

Currently, we have at least one example of a persistent supercritical accretor in our 
Galaxy - the SS433 system, where most of accreting gas is being lost in a wind. Numerical 
simulations of this system partially support the outflow-dominated scenario. However, Okuda 
(2002) states that his thorough 2D simulation fails to reproduce the outflow rate and jet colli- 
mation in SS433. In Ohsuga et al. (2005) a supercritical accretor accreting at 10 3 Eddington 
rate appears a bright (about 10 38 ergs -1 ) hard X-ray source if viewed edge-on. Ohsuga et al. 
(2005) and Heinzeller et al. (2007) calculate the structure and emergent spectra considering 
only the inner parts of the flow (R < 500Rg) not taking into account that the region considered 
is coated by accreting and outflowing matter both being optically thick. The outcoming spectra 
will strikingly differ from those calculated by Heinzeller et al. (2007). 

Optical and radio (Blundell et al. 2001) observations allow to measure the mass ejection 
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rate in the relatively slow (1000 — 2000kms _1 ) accretion disc wind seen in optical emission and 
absorption lines (Fabrika 2004) as well as the mass loss rate in the collimated mildly-relativistic 
jets launched along the disc axis. Infrared excess was used by Shklovskii (1981) and Van den 
Heuvel (1981) to estimate the mass ejection rate from SS433. No direct mass transfer rate 
measurements were ever made though it is usually supposed that it could not be much higher 
than the mass ejection rate M ej ~ lO^Moyr" 1 (Fabrika 2004). 

The details of the processes in the supercritical disc itself are unclear mainly due to the 
strong absorption in the wind. In fact even the photosphere of the wind is difficult to study 
due to high interstellar absorption (about 8 m ) making it impossible with the contemporary 
instrumentation to trace the Spectral Energy Distribution (SED) in the far ultraviolet (UV) 
spectral range were most of the radiation is expected to be emitted. The only and yet crude 
estimate of the blackbody temperature of SS433 was made by Dolan et al. (1997) and probably 
corresponds to the photosphere of the wind. Their measurements are consistent with a (2 — 
7) x 10 4 K blackbody source. 

Both observations of SS433 (Fabrika 2004) and numerical simulations (Okuda 2002; 
Ohsuga et al. 2005) support the conception that is traditionally called "supercritical funnel". 
It assumes that nearly all the matter accreted is being ejected in a form of a slow, roughly virial 
at spherization radius, dense wind. Due to centrifugal barrier two conical avoidance sectors 
with half-opening angles of 9f ~ 30° (Eggum et al. 1988) are formed along the accretion disc 
symmetry axis, filled with rarefied hot gas, which may be accelerated and collimated to form 
relativistic jets. In the inner parts of the wind this gas also can form a pseudo-photosphere 
("funnel bottom") at Ri n ~ 10 9 cm (Fabrika 2004). This radius is calculated in the consideration 
all the material ejected in the relativistic jets is uniformly distributed with respect to the polar 
angle. In case of any inhomogeneity of the jet material the inner radius becomes lower. The 
visible part of the wind may be therefore divided into three parts: funnel wall photosphere (or 
"photocone"), the outer photosphere and the inner photosphere inside the funnel. 

Another class of objects supposed to be supercritical accretors are extragalactic 
Ultraluminous X-ray sources (ULXs). Katz (1986) supposed that an object similar to SS433 
seen at low inclination angles can appear a bright X-ray source with super-Eddington apparent 
luminosity. ULXs were discovered about that time by Einstein (see Fabbiano (1989) and refer- 
ences therein). Though the nature of these sources remains unclear they are good candidates 
for the objects predicted by Katz (Poutanen et al. 2007). 

The question about how representative is SS433 among the binary systems in the su- 
percritical accretion regime in the observed Universe is difficult to answer. If one considers 
mass transfer in thermal timescale in a black hole high mass X-ray binary, the most relevant 
is the mass of the secondary that determines the timescale and hence the scaling for the mass 
accretion rate. It may be shown that in assumption the radius of the secondary scales with its 
mass as R oc M 1 / 2 and the black hole mass is close to 10 M & dimensionless thermal-timescale 
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mass transfer rate depends on the secondary mass as: 
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For the case of SS433 (M2 — 20 M & ) our estimate predicts rh ~ 2000, similar to the 
observed value 1 . If the donor star is highly evolved, accretion rate is higher and less stable 
and the phase itself is shorter. Unfortunately, initial binary mass ratio distribution for massive 
stars is poorly known but there are indications that mass ratios close to 1 are much more 
probable (Lucy 2007; Kobulnicky & Fryer 2007). The thermal timescale mass transfer rates in 
high-mass X-ray binaries are therefore likely to be highly supercritical, rh ~ 100 — 10 4 . Lower 
yet supercritical rates may take place in case of wind accretion for example in WR+BH binaries 
(Bauer & Brandt 2004). Those are likely to form a certain sub-sample of evolved ULXs with 
moderately supercritical accretion rates. 

The paper is organized as follows. In the next section we construct a simple analytical 
model describing the internal structure of an optically thick wind flow that is likely to develop 
in the case of very high mass ejection (and accretion) rate rh. Section 3 is devoted to effects 
of self-irradiation that are likely to play very important role in our funnel solution. In section 
4 we consider the emergent SEDs for arbitrary inclination taking into account self-occultation 
effects. In section 5 we test it for two sets of publicly available X-ray data. We discuss the 
implications of the model and its early testing in section 6. 

2. Structure of Supercritical Disc Winds 

In this section we study the spatial structure of the radiation density field inside the 
wind. We will use the following set of assumptions (their reliability will be discussed below in 
Section 6): (i) the flow is axisymmetric (also symmetric with respect to the disc plane) and 
stationary, (ii) the flow is optically thick to true absorption, (in) the velocity and density fields 
are not affected by energy and entropy transfer (i. e. we consider the wind already accelerated 
or accelerating/decelerating with a given power-law dependence on radius), (iv) all the motions 
are purely radial and non-relativistic, (v) the inner surface of the funnel, the funnel bottom and 
the wind photosphere are considered locally blackbody sources and (vi) for certainty we suggest 
the temperature of the bottom equal to the starting temperature of the walls. We simplified 
the picture somehow suggesting the mass is loaded in the center of symmetry. This assumption 
is violated in the inner parts of the wind. 

Effects of special relativity (namely, Doppler boosts) may be approximately accounted 
for during the calculation of outcoming spectra. We also discuss the possible influence of 

1 Acutally, there are no direct measurements of the mass transfer rate for SS433. However, the observed 
stability of the orbital period and evolution-time considerations exclude significantly higher mass transfer 
rates (Fabrika 2004; Goranskij et al. 1998). 
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relativistic effects in sections 4 and 6. 

We assume that fraction K< 1 of the accreting mass is ejected in the wind. Because mass 
ejection rate is more relevant (see below) we will use the dimensionless notation rh = M/M cr 
for the mass ejection rate where 

M cr = ~ 3 x 1(T 7 M 10 Meyr- 1 (1) 

is the critical Eddington accretion rate as intruduced by Poutanen et al. (2007). M\q here is the 
accretor mass in 1OM units. Mass accretion is characterised by raj K ~ rh, we assume K — 1 
everywhere below bearing in mind that K < 1 results in a hotter and more luminous source but 
requires higher accretion rate. The outcoming spectra depend mostly on the mass ejection rate. 
Mass accretion rate affects only the luminosity injected at the inner boundary in a logarythmic 
way. To achive a 50% change in luminosity (and about 12% change in temperature) for rh = 100 
one should assume 90% of the accreting material accreted by the black hole (K = 0.1). 

There are physical differences between the regime of highly supercritical accretion that 
we consider here and moderately supercritical accretion (rh ~ 1 — 10). At rh ~ 10 the flow 
becomes translucent and relativistic and some of our approximations are violated. 

For radial coordinate normalization we will use the "spherization radius" defined as: 
3Mk 18 „ 

R sp h = -7: — = ^-Ram. (2) 

Rq = 2GM/c 2 here is the gravitational (Schwarzschild) radius of the accretor. This value is 
proportional to the spherization radii used by Shakura & Sunyaev (1973) 



and by Poutanen et al. (2007) 

R%1 = 5/3m, (4) 

but has also non-trivial dependence on geometry. Qf = 47r cos 9 f is the solid angle of the wind. 
Normalized radial coordinates will be hereafter denoted by small letters r in contrast to capital 
R reserved for physical distances. 

We assume fixed geometry for the wind funnel, i.e. fixed half-opening angle 8f inde- 
pendent of the accretion/ejection rate. Note the difference with funnels in Polish doughnuts 
(Abramowicz 2004) where 



Of oc ^r m /r out oc e ~^ L l L ^ oc m" ' 3 , (5) 

defined by equiponetial surfaces. In contrast, when the wind is launched from a thin but 
supercritical disk at a certain radius and its velocity in the frame comoving with the disk is at 
any given starting radius proportional to the virial velocity (v = where vk is the Keplerian 
velocity) and normal to the disk surface, the half-opening angle will be a function of £ only 
(Poutanen et al. 2007): 
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Of = arctan l/J^ 2 - 1 . (6) 



Radial velocity of the wind may be uniform (v = const) in the case of large initial speed or 
virial (v oc i? -1 / 2 ) in the parabolic case when the initial velocity is close to the escape velocity 
from the wind acceleration region. We adopt a generalized self-similar scaling 



v = -y cos Of cm 1 / 2 r a (7) 

with approximately virial value at the spherization radius, a = corresponds to a constant 
velocity wind with no acceleration. In the parabolic case, when the velocity is proportional to 
the virial at any given radius, a = —0.5. 

Outside the spherization radius the gas density in the wind scales as 

and vanishes (or, at least, drops by several orders of magnitude) inside the funnel. Inside the 
spherisation radius deviations from this law are expected and mass ejection rate should depend 
on radius roughly as m e ff oc r, and hence 

nar ( - 1+o) . (9) 

We briefly analyse the consequences of this difference in density slope in section 6. 

2.1. Diffusion Equation 

The main equation governing energy transfer may be derived from the energy conserva- 
tion and Fick's laws for the thermal energy flux 

q=-DVu, (10) 

where D = c/(3np) is the diffusion coefficient and u is energy density. The first law of thermo- 
dynamics 

d (-) = Tds + ^rdn (11) 
\nj n z 

du = nTds + ^ J> ~ U dn (12) 



n 



where 



Tds = (Vq)dt (13) 

du/dt = d t u + {yV)u (14) 

dn/dt = -n(Vv) (15) 

may be rewritten in the generic form 

9 ( « + (vV)m- V(£>Vw)-7w(Vv) = (16) 

If the specific heat ratio 7 does not depend heavily on other parameters, this may be 
rewritten for an axisymmetric system in terms of enthalpy density h = as 
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\d t u + d r u - (r i+a d r u 

where 



L e = ^-d (wx6d e ) (18) 
smv 

is the polar part of the Laplace operator, and the second term in brackets accounts for adiabatic 
losses. Equation (17) is linear hence one may separate the variables. Eigenfunctions will depend 
on cos# as linear combinations of Legendre polynomials and Legendre functions of the second 
kind Q k (Jahnke & Emde 1960). 
In the stationary case 

u k = (Pk(cos6) + a k Q k (cos6))Rk(r) (19) 
where Rk is the solution of the equation 
-^e- 1 / r d r {r 4+a e 1 / r d r u)- 

_(jfc(Jb+l)-2^2l)« = (20) 

For the zeroth-order (k = 0) solution asymptotics may be derived both for large distances 
r ^> 1 where diffusion dominates, e~ x l T ~ 1 and u oc r~^ +a \ and for very low radii where 
adiabatic losses prevail and u oc r~" / ( 2+a \ Figure 1 shows the exact numerical solution for the 
radial eigenfunction Rq for a = and 7 = 4/3. It may be rather well (with accuracy better than 
2%) approximated by the following function 

F(r) = (l-e^^) 1 " (7 - 1)(2+a) r-^) (21) 

2.2. Boundary Conditions 

For a system with equatorial plane symmetry the solution will contain only even eigen- 
functions. If we choose the inner radius Ri n well below the spherization radius, equation (20) 
becomes a first-order one, and is governed by a single parameter which we assume to be the 
initial energy density. The latter may be estimated from the total central source luminosity L 
for the advective flow as 

L 



" m ' QRlv(r in )Y(r in ) (22) 
Luminosity of a supercritical disk in different models depends logarithmically on the 
dimensionless mass accretion rate (Poutanen et al. 2007; Shakura & Sunyaev 1973) as 

L = L EM (l + 0.6\n(m/K))~ 

~ 1.5 x lO^^Mio (2.8 + 0.61nm 3 ) ergs" 1 , 

where rh^ is the mass ejection rate in 10 3 M cr , X is hydrogen mass fraction in the accreting gas, 
M w is the accretor mass in 10 M units. Temperature of the flow at the spherization radius is 
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Fig. 1. Radial dependence for the zeroth-order solution for a = and 7 = 4/3 (solid line) normalized by 
the low-r asymptotics. High-r power-law asymptotics u oc r~( 2+a > is shown by a dashed curve. 

T ~ O.lm^ 3 ^ 8 keV while the gas density is n ~ 4 x 10 16 m^ 1//2 cm -3 , which leads to a rather high 
radiation to gas pressure ratio: 
aT 3 

(3 ~ — — ~ 500rh3 5/8 (24) 

For accretion rates that appear during thermal-timescale mass transfer 10 < rh < 10 5 (see also 
sections 1 and 6) radiation pressure dominates in the advective inner parts of the wind. That 
allows to equate 7 = 4/3 anywhere in the flow because the wind is radiation-pressure dominated 
in the regions where 7 is relevant. 

Energy flux is directed radially in the inner advective parts of the wind and is practically 
normal to the walls near the funnel wall surface in the outer parts. In figure 2 we show the 
two-dimesional distribution of energy density and energy flux vector field for our solution. The 
exact angular dependence of the energy influx at the inner boundary is not very important 
as it is mixed near the spherization radius, so we will restrict ourselves to the zeroth angular 
harmonic solution. 
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Funnel wall cools efficiently, therefore a simple boundary condition 

uM,) = (25) 

may be adopted. The following zeroth-order solution satisfies it and decays fast enough at 

infinity to account for the photon escape through the funnel walls 

, ( 1 — cos# l+cos0A, w , .„„. 

u = u x In — x -f Y(r) 26 

\l-cos6f l + cos# / 

uq normalization may be derived from the total luminosity (more accurately than in (22)) of 
the disk as 

U ° ~ n>Rlv(r m )Y(r m ) (2?) 

Here we define Q' = — 87rlnsin#y. Y(ri n ) ~ r i ^ 2+a ' 1 for rj n <C 1. 

2.3. Surface Effective Temperature Distribution 

Below we consider funnel walls, bottom and outer photospheres as blackbody sources 
with the temperature determined by the energy flux normal component: 

F n = a B T\ (28) 

For the funnel walls it is equal to the latitudal component: 

2L 1 

Fd = fi'sin^ / ^r^F(r m ) rl+QF(r) ' (29) 
and in the case of outer photosphere 

F r = -Dd r u = gg )L a+ 4. , r 1+a Y(r) x 



1— cosBf 1+cos 

Here r = r out should be substituted for the outer photosphere. 

Note that the fluxes are of the same order and their ratio does not depend on the radial 
coordinate. 

For the funnel walls one obtains: 

T wall = 0.038 (tt't a n 2 9 f r% a Y(r in )y 1/4 x 

x (2.8 + 0.61nm 3 ) 1/4 m3 1/2 ^ro 1/4 (^ 1+a ^(^)) 1/4 keV 
where M 10 is the mass of the accretor in 10 M units. Effective temperature scales with radius 
(for a = 0) as 

T wa ii oc r~ 5 / 12 (32) 
at lower radii and as 

T wa u oc r^ 2 (33) 

in the outer parts of the funnel. Temperature slope depends very weakly on all the parameters 
and is similar to the value p — 1/2 (if T oc r~ p ) characteristic for slim discs (Abramowicz 
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Fig. 2. Two-dimensional distribution of energy density in the optically thick wind. Thin solid 
lines correspond to logarythmically spaced (by one order of magnitude) constant energy density lev- 
els. Of = 0.4 rad ~ 23°, a = 0. Inner funnel wall is shown by a thick solid line. Arrows 
are oriented in the energy flux direction, their length proportional to the energy flux logarythm. 
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2004) and supercritical accretion discs with massive outflows, p-free multi-blackbody spectra 
are characterised by a power-law intermediate asymptotics with photon index T = 2/p — 2. 
Many ULXs have power-law spectra with r ~ 2 close to the expected in all these models. The 
similarity makes it very difficult to distinguish between the various supercritical accretion disc 
models by their X-ray spectra only. 

We assume the temperature of the funnel bottom equal to the effective temperature of 
the adjecent wall surface 

T hot = 0.038 (tt't a n 2 9 f r in )- 1/4 



7 f in) 

;2.8 + 0.61nm 3 ) 1/4 m3 1/ "M 1 " ±/ ' i keV 



Minimal inner radius Ri n may be estimated as the last stable orbit radius devided by cos Of 
as ri n = 1 /6m in R sp h units. Temperatures characteristic for the high-energy cut-offs in ULX 
spectra are about one keV and require r in > rh" 1 . It is clear that in the framework of our 
model the inner parts of the funnel should be practically transparent. If we make another 
assumption that the soft excess in ULX spectra corresponds to T ~ 0.1 — 0.2 keV temperature 
at the spherisation radius than accretion rates required are rh ~ 100. In section 5 we fit real 
X-ray spectra with our model obtaining similar results (rj n ~ 10 -3 and rh ~ 100) for two ULXs. 

Effective temperature of the photosphere may be found in a similar way. The outer 
photosphere radius may be estimated as 

l/(l+a) 



r out ~ ( ; V , I (35) 

(1 + a)Jcos9f 

in R sp h units. The photosphere is significantly non-plane-parallel, but its luminosity and spec- 
tral energy distribution may be roughly estimated as: 

(3 + a)(2 + a)( 2 +^- 1 )- 1 - 
- ~ r 2 + J (r , cos 0,L (36) 

'in 1 V in) 

T ph ~ 0.0048/ a (r*+ a Y (r in )) ^ y&^x 

x M~ 1/4 (2.8 + 0.61nm 3 ) 1/4 x (37) 

x m 3 - 3/4(1+Q) (in (^f- x ^)) V4 keV, 

•J \ \1~ cosOf 1+COS& J J 

where 

f a = (3(1 + a)) 1/2(1+Q) (2 + a )((a+«)(7-i)-i)/4 x 

x(fi/) _ (3+Q )/4(l +a )_ 

Characteristic temperatures derived above are fairly consistent with those reported by 
Poutanen et al. (2007). For a reasonably high rh = 10 3 and Of = 0.4rad the outer photosphere 
temperature is ~ 10 5 K. Together with a photospheric luminosity ~ 10 39 ergs -1 this makes 
the wind photosphere both a bright UV/optical and a bright extreme ultroviolet (EUV) source 
capable for ionizing large amounts of interstellar and circumstellar gas. We discuss the prop- 
erties of HH-regions created by supercritical accretor wind photospheres in section 6.4. Very 
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shallow temperature decline makes the spectrum flat from E ~ 0.005 keV to ~ 1 keV. The 
highest temperature value is predicted by accretion disc theory and is much higher than the 
temperature at spherisation radius (Poutanen et al. 2007): 

T max = 1.27Mr 1/4 keV. (38) 



3. Effects of Irradiation 



In the standard disc model irradiation is a minor, often negligible effect as the disc 
is thin and the fraction of emitted luminosity absorbed by the outer parts of the disc is of 
order of 0(H/R) 2 . On the contrary, in the case of supercritical funnel most of the radiation 
emitted at any point on the wall surface will be absorbed again or reflected. The probability 
of re-absorption is of the order cos#/, hence the energy balance will be significantly affected 
by irradiation. As long as true absorption dominates over electron scattering, o ^> o~t, all the 
quanta may be considered absorbed and re-emitted by the walls remaining a locally blackbody 
source. 

The energy input due to incident radiation is characterized by the flux F'(R) normal to 
the wall surface and directed inward, which may originate both from the funnel bottom and 
inner parts of the walls. In general, the incident flux may be expressed as 

X//(R')Vl(n-d)(n'.d)| 



F' 




l -dS' 



(39) 



t 6// (r) ; d* 

where the meaning of n and d vectors is illustrated in figure 3 for both wall and bottom 
irradiation. 

3.1. Problem Formulation 

In this section Cartesian and spherical coordinates are used, z axis coincides with the 
symmetry axis of the disc, funnel and jets. The remaining degree of freedom is adjusted to set 
to zero the azimuthal coordinate of the point under consideration, where the incident flux is 
calculated. This point is set by a radius vector R, 



/ sin# \ 
R = R 

\ cos6» j 

In the same way in the following subsections we define the bottom radius vector Ro and 
the radius vector R' of the variable point on the funnel walls. 



Rn — -Rn 



( o \ 


V cos6> J 



( sin # cosy? \ 
sin 9 sin ip 
cos 9 
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The geometry of funnel self-irradiation is shown in figure 3. Variable infinitesimal plate 
with area dS', radius vector R' and normal n' contributes to the incident flux in the current 
point by: 



dF'= |(n-d)(n'-d)| 



dS' 



2ixd v 

where d = R/ — R, d — |d|. Due to surface brightness invariance from distance it is convinient 
to use dulition factor w depending on geometry only and to normalize the incident flux over 
the outcoming flux -F(R) in the current point. 



w 



fo-d)(n'-d)| 

J 2 



Incident flux in the units of the outcoming integral flux (/ = F'/F) can be expressed as 
an integral over the funnel inner surface: 



/ 



'T e// (R')Vl(n-d)(n'.d) 



■dS' 



(40) 



T eff (R)J d? 

In the case of the funnel bottom we avoid integration suggesting the bottom a point-like 
source emitting as a flat plate. 

3.2. Irradiation by the Funnel Bottom 

The first heating term is simpler to handle because the source is practically point-like and 
its temperature is not affected by the funnel walls. Due to this, one can do without integration. 
Funnel bottom can be considered a spherical surface with the surface area S sp h = 2tt(1 — cos8)Rq 
or a flat circular plate of radius Rosin6 (because Rq is the starting radial coordinate along the 
wall, not along the z- axis) with the surface area S = ttRq sin 2 9. Both expressions give similar 
results when Scl. Here we assume the bottom having flat surface (see figure 3a). The photon 
source coordinate and normal unit vector are: 
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Distance vector connecting emitting and receiving points: 



d = R — Rq = Rq 



Here r = Rj Rq. Dilution factor in that case: 
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Fig. 3. Schemes explaining the geometry of irradiation effects. (a) ir- 

radiation by the bottom. (b) irradiation by the funnel walls. 



Total normalized irradiating flux from the bottom can be therefore expressed as: 
1 (r - llsin^cos 2 ^ _ A 



fi 



bottom 



(42) 



2(l + r 2 -2rcos#f 

3.3. Self- Irradiation by the Funnel Walls 

The energy input from the absorbed photons is directly calculated using equation (39). 
Surface element is an elementary area on a conical surface, alS' = 2tt sin 9 R'dR 'dip. Radius 
vectors and normals are defined. Distance vector: 



d = R' - R = R 



J 



' sin6'(l — xcos<^) ^ 

— x sin # sin <£i 
\ cos6 l (l — x 

where x = R' / R. Dilution factor in this cases takes the form: 

1 rrsin 2 0cos 2 9(1— cosip) 2 

27rR2 2((l-x) 2 +2x(l-cos ¥ 3)sin 2 6») Li 

Absorbed flux is given by an integral over x: 

''Rout/ R 



(43) 



F'(R) = 27Tsm6 / F(Rx)I{x,6)xdx 

•> Rinl R 



/R 

where I = f^ n w(x,8,ip)dtp. Integration over if is straightforward however complicated, so we 
leave it for appendix 1. Finally, normalized absorbed flux can be expressed as: 

cos 2 6> rRout/R f T eff (H') \ 4 

(44) 
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dx 




Fig. 4. Funnel wall temperature dependence on radius, normalized by r -1 / 2 . Thick solid 
line: without irradiation, thick dotted line: with irradiation effects included, three itera- 
tions. Thin solid lines represent the best power-law fits to the irradiated funnel wall temper- 
ature at larger r (> 1) and lower r (10~ 2 < r < 1), T oc r~ - 44 and T oc r ~ - 50 , correspond- 
ingly. Dot-dashed line represents the T oc r~ 3 / 4 temperature slope characteristic for standard discs. 

Irradiation is treated iteratively. Three to five iterations generally sufficient for conver- 
gence with accuracy better than 5%. As may be seen in figures 4 and 5, mostly the inner 
parts of the funnel walls (and subsequently mostly the X-ray range) are affected. Incident flux 
produced at large radii by the irradiating funnel bottom and inner parts of funnel walls changes 
as F' oc r~ 3 (see equation (42)), while the predicted temperature dependence on radius implies 
F' oc r~ 2 for the emergent flux, therefore the effect is expected. 

4. Emergent Spectra and Spectral Variability 

Here we calculate the outcoming spectra as functions of inclination and supercritical 
wind parameters. All the surfaces are considered locally blackbody sources and the observed 
spectrum (in terms of apparent isotropic luminosity vL u ) is calculated by integrating over all the 
visible surfaces applying visibility conditions. The most important effect is that at inclinations 
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Fig. 5. Irradiation effects on the outcoming spectrum. Solid lines show the spectral energy distribution 
of a face-on supercritical funnel for different number of iterations (0, 1, 2 and 3). Funnel walls without 
irradiation give the spectrum shown by dashed lines. Dotted lines correspond to the contributions from 
the bottom (at high energies) and the outer photosphere. Parameters are the same as in figure 2. 

i^jOf funnel bottom and some part of the funnel wall become obscured from the observer. 

Doppler boosts are accounted for approximately. We consider funnel walls moving at 
an angle 8f with the line of sight (that is exactly true for a face-on object), and the material 
consisting the funnel bottom moving at an angle i. The velocity is determined according to 
equation (7) and is considered equal for the funnel bottom and the wind at r in . Photon energies 
are modified by a factor: 

1 



7(l-/3cos0) 



(45) 



where f3 = v jc is the dimensionless velocity, 7 = 1/y 1 — (v/c) 2 is the Lorentz factor, and 9 is Of 
for the walls and % for the bottom. Photon numbers are modified by a factor 5, and the time 
compression factor is not relevant here because the flow is stationary. Apparent luminosity 
L E changes proportionally to L E (E) oc 5 2 L e j \e/5), where L E is the apparent luminosity 
calculated without relativistic effects. Note that Doppler boosts become important only when 
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the velocities in the inner parts of the flow are relativistic (either when fn < 10 or for a ^ —0.5). 
See also discussion in section 6.2. 

In figures 4 and 5 we compare the temperature profiles and SEDs with and without 
irradiation. We adopt a = 0, rh = 10 2 , Ti n = l/6m = 1.7 x 10 -3 , Of = 23°. The temperature 
dependence on radius is close to a broken power-law in both cases, and the net effect of irradi- 
ation is in altering the mean temperature in the innermost parts of the funnel by about 20%. 
In figure 5 SEDs are calculated for zero inclination and different number of iterations used to 
account for irradiation effects. Evidently, only the funnel walls is affected by irradiation that 
alters the flux by a factor of 2 — 3. Spectrum is practically flat in the EUV/soft X-ray region 
but curves near the T sp h ~ 0.1 keV. High-energy cut-off is present at several keV if the inner 
radius is rj„ ~ 1/rh. 

4-1- Self-Occultation 

Let i be the inclination of the funnel (angle between the symmetry axis and the line 
of sight). For any given radial coordinate at the funnel wall surface three cases are possible: 
the annulus is fully visible (this is possible only for i < Of), invisible or partly visible. If the 
inclination is larger than Of for different radii either second or third case takes place for every 
annulus. Visibility may be quantitatively described by a factor y defined as the visible part of 
the given annulus (having constant distance r from the center). It may be determined by an 
analytical ray-tracing method. 



Here 1 is a unit vector directed towards the observer (we use the same scheme as in the previous 
section), R is the radius vector of the intersection point between the funnel surface and the ray 
starting from r and directed towards the observer. Boundary case when the intersection point 
has the radial coordinate equal to the photosphere radius is of interest here. Normalizing over 
the outer radius (x = r/r out ), one can express the vectors as follows. 



r = R + sl 



/ smOfCosA 



r = x sinOfsinA 



\ cos9f 
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It is easy to express two of the unknowns (s and Aq) from the others and obtain the 
solution for A (azimuthal coordinate of the starting point): 



If the annulus is partially visible y may be calculated as the length of the arc (devided 
by 2tt) connecting the two points defined by the solutions of the equation above. Finally: 



4-2. Implications of Self- Occupation Effects 

Watarai et al. (2005) studied self-occultation effects for fat accretion discs coming to 
an evident conclusion that the spectrum becomes softer at large inclinations. Heinzeller et 
al. (2007) comes to a similar conclusion considering the results of 2D radiative hydrodynamic 
simulations by Ohsuga et al. (2005). In our case this effect is even more profound (see figure 6): 
at inclinations i ~ 9f the X-ray component abruptly disappears (because only the inner parts 
of the funnel contribute significantly to the X-ray range) and only the EUV component may 
be observed. 

Due to that reason some supercritical accretors (those viewed at large inclinations) will 
not show the ULX phenomenon (see figure 7) but remain luminous UV sources. For 8f = 0.4rad 
the number of orphan ultraviolet sources of similar nature is about cos#// (1 — cos#/) ~ 10 times 
higher than the number of ULXs. 

Observations of these sources in the EUV are complicated by neutral gas absorption. In 
Abolmasov et al. (2008) we discuss the possibility of far-ultraviolet observations with GALEX 
(Martin et al. 2005) and HST in the ultraviolet. Supercritical accretors prove to be difficult 
but still reachable targets for GALEX and HST, appropriate for pointing observations. 

We do not know yet how much there is similarity between ULXs and SS433. SS433 
exhibits jet/disc axis direction variations with a super-orbital precession period about 160 d . 
Similar effects are observed for some other X-ray binaries (Clarkson et al. 2003). There are 
indications that the supercritical disc (and consequently the wind with the funnel) follows the 
motions of the jets. 

A nearly face-on supercritical accretion disc may show periodical X-ray variability by 
one-two orders of magnitude when 8f is less than the maximal and higher than the minimal 
value of i. Precessional variability should have a distinguished light curve with a flat maximum 





(46) 



where: 
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Fig. 6. Apparent isotropic luminosity as a function of inclination. Solid line represents bolo- 
metric luminosity. Dashed, dot-dashed and dotted lines represent harder X-rays (0.4—20 keV), 
softer X-rays (0.1 - 0.4 keV) and EUV (0.01 - 0.1 keV). Vertical lines mark f and tt/2 - 9 f . 

and continuous flux change near the minimum light. In figure 8 we present the predicted X- 
ray spectra as functions of precessional phase for two sets of parameters. We assume here a 
simplified version of the kinematical model applied to SS433 (Abell & Margon 1979): 

cosi = cosi^cosio + smiAsmi cos(27rip), 

where io is the inclination with respect to the precession axis, %a is the amplitude of inclination 
variation, and is the phase of super-orbital period. 

The only ULX exhibiting variability that may be considered super-orbital is X41.4+60 
in M82 (Kaaret et al. 2006). RXTE X-ray flux varies by ~ 50% for this source with a 62 d 
period. Irregular variability on similar timescales is much more ordinary among ULXs (La 
Parola et al. 2001) possibly indicating that the funnel shape itself changes rapidly enough to 
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Fig. 7. Inclination effects on the outcoming spectra. Apparent values of vL v are shown 
for i = 0, 25, and 90 degrees as a solid, dashed and dotted lines, respectively. 

level the effects of precessional variability. 
5. Comparison with Observations 

5.1. Data on ULXs 

XMM-Newton datasets on two X-ray sources, NGC4559 X-7 and NGC6946 ULX-1 were 
analysed. We used archival XMM-Newton EPIC (Turner et al. 2001) data (MOSl,2 and PN 
detectors) obtained on 8 October 2006 (observation ID 0152170501) and 25 June 2004 (observa- 
tion ID 0200670401), correspondingly. All the data were reduced using standard XMM-Newton 
Science Analysis System (SAS) procedures. Response files were made using SAS tools rmf gen 
and arf gen. We set flag= for PN data. We used pattern < 4 for PN and < 12 for MOS data. 

NGC4559 X-7 is known as a "supersoft" bright ULX (Cropper et al. 2004). It does not 
coincide with a bright stellar or nebular optical counterpart. However, there are indications 
that the source is connected with ~ lOMyr old stellar population (Soria et al. 2005). 

NGC6946 ULX1 (Roberts & Colbert 2003) is known to be a source of moderate luminos- 
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Fig. 8. Spectral variability of processing supercritical accretion disc funnel. Lines of con- 
stant flux (EFe) as a function of energy and precessional phase are given. Flux is normal- 
ized over its maximal value. Of = 22.9 deg, rh = 10 3 , a = 0. Upper diagram corresponds 
to io = 20 deg and ia = 20 deg, lower has io = 45 deg and ia = 20 deg. In the latter case 
funnel bottom is never visible and funnel walls always remain at least partially occulted. 
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ity Lx ~ 3 x 10 ergs - with an X-ray spectrum usual for ULXs. Stellar optical counterpart was 
detected by HST'va the visible range (Blair et al. 2001) as a relatively bright star (My ^ — 7 m ). 
As many other ULXs, NGC6946 ULX1 coincides with a bright optical nebula. Abolmasov et 
al. (2008) analysing the emission spectrum of the nebula conclude that the central source may 
be a hard EUV source with temperature T ~ 10 5 K and luminosity L ~ 10 40 ergs -1 , roughly 
consistent with what one may expect from a m ~ 10 2 — 10 3 supercritical accretor. 

Three models were used to fit the data: standard disc (Mitsuda et al. 1984) + power 
law, p-free disc (Mineshige et al. 1994) and self-irradiated multi-color funnel model presented 
here. Efficiency of the p-free model for ULXs was already shown by Vierdayanti et al. (2006) 
and other works. Because the temperature of the funnel walls decays in a practically power-law 
manner (with p ~ 0.5, see figure 4) one should expect p-free model to produce SEDs close to 
the SEDs predicted by our model. However, the temperature distribution required is very flat, 
so we have shifted the lowest allowed p from 0.5 to 0.05. For each object all the three extracted 
spectra were fitted simultaneously. Spectral ranges 0.1 — 12keV for PN and 0.1 — lOkeV for 
MOS data were used. Photoelectric absorption by a solar metallicity material was included as 
a free parameter in all the models. 

For the supercritical funnel model used for X-ray data fitting we fixed the velocity 
exponent (a = —0.5), the outer radius {r out = 100) and mass ejection rate (to rh = 10 4 ) and the 
inclination to 0°. For fixed T in , r in and r out values the only effect of the mass accretion rate is 
in changing the gas velocities. As long as v{r in ) -C c, the mass ejection/accretion rate m does 
not affect the shape of the X-ray spectrum. 

Low values of a ~ 0.. — 0.5 provide better fits to the data. However, we fix the parameter 
to avoid degeneracy with Tj n and ri n . a = —0.5 is expected if the velocity is proportional to the 
virial velocity at any given radius. In ballistic approximation (when the particles of the wind 
are first rapidly accelerated and then move along hyperbolic trajectories in the gravitational 
field of the accretor) a ~ —0.5 in the inner parts of the flow and approaches at larger radii. 

In table 1 we present the results of spectral fitting. It may be seen that in both cases a 
two-component standard disc + power law model gives the best results because of being capable 
to fit the harder part of the spectrum that is difficult to handle using only thermal models with 
exponential high-energy cut-offs. 

Because the best-fit temperature in both cases is rather high Tj„ ~ 1 — 2keV, we expect 
the inner parts of the funnel to be practically transparent and Ri n ti be close to the last stable 
orbit radius. In this assumption, Tj n may be directly converted to the mass accretion rate. If 
we equate the inner radius to the last stable orbit for a Schwarzschild black hole, 

cos 9 f R sph 1 
18 R G Qr in 

The estimated dimensionless mass ejection rates are therefore about 100 for NGC4559 X- 
7 and about 300 for NGC6946 ULX-1, correspondingly. The latter value is consistent by the 
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order of magnitude with the mass accretion rate estimate for NGC6946 ULX-1 resulting from 
the optical spectroscopy of the nebula MF16 associated with this X-ray source (Abolmasov et 
al. 2008). 

The well-known high-energy curvature (Dewangan et al. 2005) can be explained by 
Ti n ~ 1 — 2.5keV. If a hardening factor ~ 1 — 2 is used (see also section 6), the inner temperature 
may be similar to the expected maximal temperature for a critical accretion flow around a 
~1OM black hole given by equation (38) or even slightly lower. High inner temperatures like 
lOkeV also result in acceptable fits. Model is insensitive to r out save for the cases when this 
parameter is < 1 therefore we fixed the parameter to 100. 

Low values of rj n appear to be a real feature of ULXs. Funnel interior is expected to be 
practically transparent to the X-rays down to several gravitational radii, having rather deep- 
lying bottom pseudo-photosphere (or none at all). Observations in a broader spectral range 
are needed to distinguish between the thermal radiation from the immediate vicinity of the 
compact object and comptonisation effects. 

Though the model fits the data in a quite acceptable way its parameters are poorly 
constraint. The spectral shape does not depend very much on the actual accretion rate and 
velocity law. There are numerous ways of making the spectrum harder: taking into account 
the difference between mass accretion and ejection rates, applying comptonisation effects etc. 

6. Discussion 

6.1. Limitations of the Model 

Our calculations do not account for radiation feedback on the dynamics. That is not 
only wind acceleration within the regions of interest but also evaporation from the funnel walls 
and development of instabilities in the outer parts of the funnel in the strong radiation field 
of the inner parts. The instability may resemble that of accretion discs with irradiation (see 
Mineshige (1993) and references therein) but the wind is unable to influence the irradiating 
source and its role in developing the instability is purely passive. 

The structure of the inner parts of the flow may be much more complicated than we 
assumed. At r ~ 1 the source of the wind becomes spatially resolved and our radial approxima- 
tion fails. Besides this, the angular distribution of the energy influx may be more complicated 
and higher angular harmonics may appear. 

Qualitatively the effects of non-zero size mass loading region may be considered as fol- 
lows: let us assume that mar. The main radial equation for k = takes the form: 



It is easy to check that the equation allows two power- law solutions (u oc r a , where a = 



—0.5 ( (1 + a) ± J(l + a) 2 + 47(2 + a))) corresponding to inward and outward diffusion. For 
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r< 1 and a = —0.5 energy density decreases as u oc r~ im instead of u oc r -2 . The difference 
increases for higher a. 

We do not account for the emission of the hot interior of the funnel as well as for thermal 
comptonisation effects in the wind and pure reflection. All these effects are likely to harden 
the outcoming spectra making inner temperature estimates from the observational data shifted 
towards higher values. Hardening factor Th — 2.6 measured for standard discs (Borozdin et al. 
1999) together with relativistic effects may lead to the inner temperatures about 1 keV appear 
as several keV. 

Reflection and absorption by moving partially ionized gas may simulate the soft-excess 
observed in many ULX spectra (Stobbart et al. 2006) in a way it was proposed by Gierlihski 
& Done (2004) for AGNi and recently by Gongalves & Soria (2006) for ULXs. Understanding 
the structure of the outer photosphere of the wind requires more complicated modelling taking 
into account both significant non-sphericity of the outflow and various opacity sources. 

6.2. Relativistic Effects 

Broadly speaking, relativistic effects are expected to change photon energies and fluxes 
by factors close to 5 = (7(1 — /3cos#/)) -1 . Here, f3 and 7 correspond to the local dimensionless 
velocity and Lorentz-factor. The temperature of the X-ray spectrum is additionally increased 
by: 

y = 5-l~5x 10~ 3 cos 3/2 6 f m^ 1/2 r° (47) 

The effect becomes significant if the velocity is relativistic in the inner parts of the flow 
(that corresponds to a < —0.5 in the above formula). The total observed luminosity is altered 
by a factor <5 3 (S 2 factor from geometrical reasons and one S from the increased energy of the 
photons). For rh = 10 3 the effect is about one percent if the velocity at the spherisation radius 
is considered. However, if the inner parts of the wind have mildly relativistic velocities, the 
X-ray part of the spectrum appears about 2 — 3 times brighter for (3 ~ 0.3, 6f ~ 20°. 

The funnel bottom is the most likely part of the flow to be affected by relativistic effects. 
The gas inside the funnel cone is supposed to move with mildly relativistic velocities (Eggum 
et al. 1988), approximately towards the observer at low inclination angles. The observed colour 
temperature changes by a factor of <5~ l+/3~ 1.1 — 1.5. The total observed luminosity becomes 
several times higher. Our model allows to include the relativistic Doppler effect in calculations 
accurately for i = (see section 4). At non-zero inclinations the deviations from our approxi- 
mation (that all the funnel wall material moves at 6f angle with respect to the line of sight) 
are of the order 0(/3smi) 2 . 

Relativistic aberration is also able to dump irradiation effects. The X-ray part of the 
spectrum will be the most affected. Irradiation becomes considerably smaller if the radiation 
emitted by the moving gas becomes anisotropic enough. However, for mildly relativistic flows 
(P < 0.5) the effect is not so severe compared to the effects mentioned above. For example, 
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the irradiating flux from the funnel bottom changes for r ~ r ~ 1/Qrh and a = —0.5 by about 
7 3 — 1 ~ 20 — 30%. At large radii the effect is smaller and may change sign at certain radii. 

6. 3. Observational Predictions for ULXs 

There are several effects expected if our estimates are correct: (i) supercritical accretors 
viewed at low inclinations (i < Of) will be seen as ULXs (L > 10 39 ergs -1 ); (ii) independently 
of its inclination a supercritical accretor is a luminous UV source; (Hi) predicted X-ray spectra 
of supecritical accretors viewed face-on are similar to p-free model spectra with p ~ 0.4 — 0.6; 
(iv) because of high EUV luminosities supercritical accretors should ionize the wind above the 
photosphere and establish Stromgren zones; (v) if the accretion disc precession characterising 
SS433 is usual for supercritical accretors viewed nearly face-on, strong X-ray flux modulation 
is expected for at least a number of sources with "super-orbital" periods like tens and hundreds 
of days. 

Most of the observational properties of ULXs are naturally explained in our model. There 
is observational evidence that at least some ULX nebulae are powered by photoionisation from 
the central source having luminosity comparable with the apparent X-ray luminosity. IMBH 
binaries should have difficulties in ionizing the surrounding gas unless the mass of the IMBH 
is very high, > 10 4 M Q (see discusion in Abolmasov et al. (2008)). 

For very high accretion rates rh ^ 10 3 the outer photosphere is very large, comparable 
to the probable size of the binary system. Applying a = and Of = 0.4 rad results in physical 
radius values: 

R out ~ 4 x 10 13 m3 /2 cm. (48) 

However in a real high-mass binary conditions the structure of the outflow is preturbed 
by tidal forces and the outflowing gas itself rapidly recombines (see next section). Because the 
flow is essentially supersonic and perturbed by a strong non-axisymmetric potential it is likely 
to become highly inhomogeneous. The equatorial outflow of SS433 is an example of such kind. 
The photosphere size in the optical then saturates at a radius of the order 10 12 cm preventing 
wind photospheres from becoming "red hypergiants" with very high infrared luminosities. The 
effect may become significant starting from fa ~ 100. 

Further understanding of supercritical accretor winds will require methods used for stel- 
lar atmosphere calculations. Rosseland mean for n ^ 10 10 cm~ 3 and T ~ 10 4 — 10 5 K is very 
close to k t (Iglesias & Rogers 1996) but in certain spectral ranges the wind should be less 
transparent. An edge-on ULX will mimick an OB-hypergiant with dense and fast wind or a 
low-temperature hydrogen-rich WR. 

6.4- Photoionized Nebulae 

In the case of highly supercritical accretion UV and EUV spectra may give much more 
information about the mass ejection rate than the X-ray properties. In figure 9 we show the 
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dependence of H and He + -ionising fluxes (1 — 1.5 and 4 — 20 Ry ranges, respectively) and 
corresponding luminosities of recombination emission lines on mass accretion/ejection rate. 
Here we suggest that the HII regions have fixed temperature T= 10 4 K and density n = 100cm~ 3 . 
Atomic data were taken from Osterbrock & Ferland (2006). One may see in the figure that the 
high emission-line luminosities observed in ULX nebulae may be well explained by high EUV 
luminosities of the central sources. The HellA 4686 / H/3 ratios predicted by the model are 
close to the high HeIIA4686 / H/3 ratios ~ 0.2 measured for some high-excitation ULX nebulae 
or the inner high-excitation parts observed in some of the ULX shells (Abolmasov et al. 2007). 

The situation becomes more complicated if one takes into account absorption in the 
wind optically thin to electron scattering. Density at the wind electron-scattering photosphere 
is 

n = — -= — ~ 3 x 10 9 m7 3/2 cm- 3 (49) 
VfRlutV 

Recombination will occur in a layer having thickness: 

an (50) 

~ 2 x 10 11 m 3 M 10 r 4 " 1/2 cm, V ; 

where a ~ 2.6 x 10~ 13 T 4 _1 ^ 2 cm 3 s~ 1 is the effective recombination rate and T 4 is gas temperature 
in 10 4 K units. The size of the recombination region is close to or smaller than the size of 
the outer photosphere. Actually that means that the outer photosphere of the wind will be 
optically thick to the EUV radiation of the central source. 

Hard radiation from the pseudo-photosphere may however ionize the gas. Because the 
density of the wind falls off rapidly two regimes appear: ionized and neutral wind. First case 
makes the photosphere of the wind a UV object with high bolometric correction and a compact 
HH-region coincident with the X-ray source. In the latter the photosphere is much cooler and 
mimicks a hypergiant with broad emission lines. 

The ability of the central source to ionize the wind may be calculated as follows. 
Recombination rate integrated over the outer wind is 

f Rmax 

1=1 aQ f n 2 (R)R 2 dR (51) 

J Rout 

Rmax here is the radius of Stromgren zone ionized by the central source. / may be used as an 
estimate for the number of hydrogen-ionizing quanta intercepted by the wind, if R max is set to 
infinity: 

S ~ 2 aR 2 ° m 3 / 2 ~ 1.2 x 10 47 mf s" 1 (52) 
cos^cosfc*/ a j. 

This recombination rate is usually lower than the quanta production rate, affecting only the 
highest accretion rates. In figure 9 we present the ionizing quanta production rates for different 
mass ejection rates. ULX nebulae may appear even brighter due to two additional effects: 
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Fig. 9. Numbers of H (solid curve) and He + (dotted) -ionizing photons produced by our 
percritical wind model (upper panel) and the luminosities of relevant recombination lines 
Ha and H/3 (solid) and HeIIA4686 (dotted). a = 0, r in = l/6m and 8 f = 23° is 
sumcd. Thin straight solid line shows the number of quanta absorbed in the wind. 
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"stripping" of the outer wind photosphere (the effect proposed in previous section) and higher 
accretion power due to K < 1 . 

7. Conclusions 

Optically thick wind with irradiation is capable to explain the SEDs of ULXs in the 
standard X-ray band and even the high-energy curvature that is difficult to explain in the 
framework of unsaturated comptonisation cool-disc IMBH model (Dewangan et al. 2005). High- 
energy cut-off is predicted to appear at several keV. Higher observed values of T in ~ 1 — 2keV 
may be explained by applying a hardening factor ~ 1 — 3 similar to those predicted for accretion 
discs in X-ray binaries. 

Outcoming spectra observed at low inclinations are similar to the spectra of slim discs 
and resemble p-free model spectra with the p parameter close to 0.5. X-ray spectra of known 
ULXs are approximated equally good by our model, p-free and standard disc + power law two 
component models. However the parameters are poorly constraint supporting the idea that 
the properties of the X-ray spectrum depend rather weakly on the accretion disc and wind 
parameters. Relatively high inner temperatures argue for the funnel interior to be transparent. 
In this assumption, r in parameter may be used to estimate the mass ejection rate that appears 
to be of the order m ~ 100 — 300 for the two sources analysed. 

We stress the extreme importance of irradiation effects providing mild geometrical col- 
limation of the observed X-ray radiation. In a simple assumption of local absorption and 
re-radiation of the absorbed energy we show that the temperature of the funnel wall surface 
is altered by about 20% in the inner parts of the funnel, and the outcoming apparent X-ray 
luminosity becomes about 2 — 3 times higher. 

Photoionized nebulae are likely to be formed around supercritical accretors. Ionizing 
quanta production rates suggest that in most cases a supercritical accretor is capable to produce 
a photoionized HH-region with bright optical emission line luminosities ~ 10 37 ergs -1 but higher 
luminosities may appear as well. 

The work was partially supported by the RFBR/JSPS grant 05-02-19710. 
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Table 1. Best fitting results for the two selected ULXs. Errors correspond to 90% confidence range. 



NGC4559 X-7 NGC6946 ULX-1 



standard disc + power law 



Nh, 10 22 cm- 2 


219° S 3 nn 
—0.02 


44+R-Sf 

"■^^—0.05 


T m ,keV 


n irrH-0.010 

u.ioo_ 009 


U.lOD_ 019 


standard disc normalization 




fi9n +2000 
DZU -200 


r 


o 9/1+0.05 
z - z ^-0.05 


o 4fi +o.i5 


power law normalization 


2.3i°,; 14 x lO" 4 


3.31° i x 10- 4 


X 2 /DOF 


648/673 


523/504 


L modeh 10 39 ergs^ 1 


9.6jl?; 4 


o i+3 
a - i -1.0 



p-free disc 

N H , 10 22 cm- 2 0.0407t°;°°S£ 0.293jl° : °? 6 
T m ,keV 0.64;°; 22 2.3l°; 4 



v 460+ 004 4 04+ - 01 

P u -^ uu -0.005 u -^ u ^-0.01 



-0.005 

p-free disc normalization 1.4 jig J x 10~ 4 3.3±§_ 4 x 10~ 4 



x 2 /DOF 734/674 665/509 
j e/ ,10 39 ergs- 1 9jJ 3±1 



self-irradiated multi-color funnel 



N H , 10 22 cm- 2 


n lfi7 +0.026 
u - 10 ' -0.015 


9Q^-0 ° 16 


T in ,keV 


1 41+0.16 
i '^ i -0.12 


1 q+0.3 


' in 


1.6+?;? x 10- 3 


6±| x 10- 4 


0,deg 


12 2 +0 - 7 

iz - z -0.6 


19±l 


normalization 


30i« 2 


74+J , 


X 2 /DOF 


660/673 


604/504 


L m odeh 10 39 ergs" 1 


10.5±2 


3t 2 
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Appendix 1. Integration over ip 



The integral over (p can be calculated as follows: 



/7T PTT (~\ COS (Dl 
w(R'/R,8 f ,<p)d<p = —xsm 2 9 f cos 2 9 f / ■) 
-7T 27r j-tt (a + frcosy?) 



Where a = 1 — 2xcos 2 6>/ + x 2 , b = —2xsm 2 8f 



i = ^'8,^0, (/.-. t-^^ - 2/:. (^=5^+ /_", (^S^ # 
= W/ ((i - S) ^^j. ^ - 2 (i + f) 



The two integrals can be expressed as follows (for details see for example Dwight 1961 
or any other table of integrals): 

^ (a+bLo f y d<P = a+bLe f d( P = ( a *-£?/* 

rn cosy j _ b r^r 1 7 _ 2ivb 

J — (a+bcose f y ^ ~ « 2 -* J-* a+bcose f ~ (a 2_ fe2) 3/2 

Finally, the integral value becomes: 

T _ xsin 2 9 f cos 2 e f f 1 _ / q+fc g-2b 



b 2 

f2. 



*(1 

277X00120/ i 7 1 |l-x| l-2x(l-3sin 2 6» / )+x 5 

v /l-2zcos(26» / )+x 2 l-2xcos(20 / )+a; 2 
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